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By comparing neutrino fluxes and central temperatures calculated from 1000 
detailed numerical solar models, we derive improved scaling laws which show 



■ how each of the neutrino fluxes depends upon the central temperature (flux 

(N 



oc T"^); we also estimate uncertainties for the temperature exponents. With 
the aid of a one-zone model of the sun, we derive expressions for the tem- 
perature exponents of the neutrino fluxes. For the most important neutrino 
fluxes, the exponents calculated with the one-zone model agree to within 20% 
or better with the exponents extracted from the detailed numerical models. 



^ . The one-zone model provides a physical understanding of the temperature 

dependence of the neutrino fluxes. For the pp neutrino flux, the one-zone 
model explains the (initially-surprising) dependence of the flux upon a nega- 
tive power of the temperature and suggests a new functional dependence. This 
new function makes explicit the strong anti-correlation between the '^Be and 
pp neutrino fluxes. The one-zone model also predicts successfully the average 
linear relations between neutrino fluxes, but cannot predict the appreciable 
scatter in a Acpi/cpi versus A(f)j/(j)j diagram. 
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I. INTRODUCTION 



Deciphering the solar neutrino problem offers the combined challenge of understanding 
the structure of the solar interior and understanding the nature of neutrino interactions. 
The consensus view at present, in part based upon temperature scalings discussed in this 
paper, is that the measured solar neutrino fluxes reported in the four operating experiments 
cannot be explained by hypothesizing changes in standard solar models (SSMs). The most 
plausible explanations, with the currently available data, require some extension of the 
standard model of electroweak interactions [T|-TB|. 



The long-standing discrepancy between the observed and the predicted neutrino fluxes 
has motivated the study of many non-standard solar models, which are in most cases ad hoc 
perturbations of the standard solar model. For many of the proposed changes of SSM input 
parameters (e.g. nuclear cross sections, element abundances, and opacities), the predicted 
neutrino fluxes are approximately characterized by a single derived model parameter, the 
central temperature, T. For small variations of input parameters, the neutrino fluxes and 
the central temperature of a detailed solar model can be related by a power law of the 



form 1 7 



(f)(xT , rn = — — -. (1) 
rflogT 

For the fundamental pp neutrinos, m ~ — 1 and for the important neutrinos, m ~ +20. 

The temperature dependences indicated in Eq. (D) are obtained from precise calculations 
with complex stellar evolution codes that solve coupled partial differential equations. The 
results of these calculations are well known and have been used in many previous analyses of 
the implications of solar neutrino experiments. Our goal is to show how these simple depen- 
dences result from the basic physics of the problem and to what extent the parameterization 
in terms of a central temperature is sufficient to characterize the neutrino fluxes. The most 
initially suprising result of the stellar evolution calculations is that the magnitude of the pp 
neutrino flux dependence depends inversely upon the value of the central temperature. We 
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shall see that even this result has a simple, quantitative explanation in terms of the nuclear 
physics of the energy generation process. 

The scaling with the central solar temperature can be used to evaluate neutrino fluxes 
for small deviations from the standard solar model. If a model is known to have a slightly 
different central temperature than the SSM, the neutrino fluxes can be estimated without 
detailed numerical calculations. Many classes of non-standard models (involving, e.g., rapid 
rotation in the solar interior, some variations in nuclear cross sections, or the existence of 
a strong magnetic field in the solar core), reduce the central temperature. A quantitative 
determination of the reduction in the central temperature implied by a specified change in 
the input physics requires detailed numerical modeling. With the aid of the temperature 
scaling laws, neutrino fluxes from non-standard solar models can be investigated over large 
parameter ranges. 

Using previously determined scaling laws of neutrino fluxes with central temperature, T 
jr7|JT8[] , several authors P,[7HlOi have studied non-standard solar models and have compared 



the predicted neutrino fluxes with the available experimental data from the four operating 
solar neutrino experiments, Homestake, Kamiokande, GALLEX, and SAGE |TP|-^. These 



authors [p|,[7|-pil[| show that it is impossible to reconcile the data from the four operating solar 
neutrino experiments with the neutrino fluxes predicted by changes in T. They conclude 
that non-standard solar models which have different central temperatures than the standard 
model are unlikely to solve the solar neutrino problem. 

In the applications described above, the conclusions depend to some extent upon the 
extrapolation of the temperature scalings to a larger temperature range than was covered in 



the original numerical models from which the scaling laws were derived |17]. We note that 



Castellani, Degl'Innocenti, and Fiorentini ^ (see also references P,p!0||) find good agreement 
between the neutrino fiuxes calculated from their non-standard, but detailed, solar models 
and the fiuxes obtained by scaling with respect to the central solar temperature. 

The analysis in the present paper provides a physical justification for the use of the 
temperature scalings, even for relatively large changes in the central solar temperature. 



We present improved determinations for the temperature exponents and estimates for their 
uncertainties. We also give exponents for four minor neutrino fluxes for which temperature 
dependences were not previously available. By construction, the derived scaling law for the 
fundamental pp neutrino flux is consistent with the observed solar luminosity. 

Our results are complementary to the powerful numerical techniques of Castellani et 
al. |l7|,|9|JlO|l , who have shown that several non-standard solar models have a homologous 



temperature dependence, and to the insightful physical analysis of Bludman who has 
argued that it is a reasonable approximation to regard the solar interior as a single region 
that can be largely described by a single parameter, the central solar temperature. Taken 
together, these previous investigations provide strong motivation for taking seriously the 
predictions of a single-zone solar model. 

In what follows, we shall refer frequently to different fusion reactions in the pp chain. 
For convenience, we summarize in Table | the principal reactions in the pp chain. 

In Sec. II of this paper, we report scaling laws (with uncertainty estimates for the expo- 
nents) for all the solar neutrino fluxes that were included in the 1000 numerical solar models 



calculated in the Monte Carlo study of Bahcall and Ulrich [|T7|. We present, in Sec. Ill, a 
simple one-zone model for the sun which accounts well for the numerically-derived scaling 
laws. This model motivates the use of the new functional form for the temperature depen- 
dence of the pp neutrino flux. We discuss in Sec. IV the Fogli-Lisi sum rule on the 
temperature exponents. In Sec. 0, we display the correlations found between the values of 
the principal neutrino fluxes in the 1000 solar models and show to what extent the one- zone 



model can account for these correlations. Finally, in Sec. ^ we summarize and discuss our 
main conclusions. 

II. TEMPERATURE EXPONENTS: DETAILED SOLAR MODELS 

We have derived power-law scaling relations with the central solar (model) temperature, 
T, for the neutrino fluxes from the 1000 numerical solar models of reference [^. The 
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numerical models are calculated with a stellar evolution code (the same code each time) 
that typically uses a few hundred mass zones. For these precise solar models, the principal 
input parameters, (nuclear cross sections and element abundances) were sampled within 
their ranges of uncertainty. The average dependences upon the central temperature T of the 
neutrino fluxes calculated for these models are represented with reasonable accuracy |T^ by 
power-law relations. 

Some particle physicists have objected to the term '1000 standard solar models' to de- 
scribe this collection of numerical solar models on the grounds that the same underlying 
theoretical model (stellar evolution theory and the standard electroweak model) is used to 
calculate all of the solar models. For their comfort, we have avoided in this paper describing 
the set of models as '1000 standard solar models' and instead refer to them as '1000 detailed 
solar models' or '1000 numerical solar models.' 

In Figs, m and ^ and in Table |I|, we show temperature exponents that were derived 
by minimizing the residuals in power-law fits of the fiuxes versus T. The power-law repre- 
sentations generally describe well the dependences of the neutrino fiuxes upon the central 
temperature of the solar model. 

The temperature exponents obtained here by minimizing the residuals for the power law 



fits are in approximate agreement with the earlier values |117| , [r8| obtained by best visual fits 
( "chi-by-eye" ) of power-law temperature dependences to the distribution of neutrino fiuxes. 
The new (previous) temperature exponents of the neutrino fiuxes are: m{pp) = — 1.1(— 1.2); 
m(^Be) = 10 (8); and m{^B) = 24 (18). 

For the first time, scaling exponents derived from the Monte Carlo experiment are given 
here for the less-numerous neutrinos from pep, ^^N, ^^O, and ^''F. In previous applications 
of the scaling laws, it was necessary to guess the temperature exponents of these four minor 
neutrino fiuxes based upon analogies with the published exponents for the three dominant 
neutrino fiuxes. 

We estimate uncertainties in the values of the scaling exponents by determining the 
indices twice, once by minimizing the residuals in the best fit to the neutrino fiuxes and once 



by minimizing the residuals in the best fit to the central solar temperature. This method 
attributes all of the scatter in the fits to either the flux or to the central temperature; 
the calculated range provides a reasonable estimate of the plausible range of the scaling 
exponents. 

Table |T| presents the best-estimates for the scaling exponents and their estimated errors.^ 
We also show in Table |l| the semi-analytic results obtained in the following section using 
the one-zone solar model. We choose to represent the pp flux with a new functional form, 
(ppp oc [1 — 0.08(T/Tssm)'"], rather than (as has become conventional) (j)pp oc T™, for the 
reasons described in the next section. 

III. TEMPERATURE EXPONENTS: ONE-ZONE MODEL 

Using a static one-zone model of the present-day sun, we derive in this section approxi- 
mate scaling laws that agree satisfactorily with the results obtained by detailed evolutionary 
calculations of numerical solutions of multi-zone solar models. We assume a fixed tempera- 
ture and matter density for the one- zone solar model. 

With this extremely simplified model, we do not expect to obtain accurate values for the 
temperature scaling laws. However, we shall see that the temperature exponents that are 
obtained agree surprisingly well with the scaling laws derived from the precise evolutionary 
solar models. 



^The residuals were minimized, as presented in Figs. || and ||, in logarithmic space. For example, 
the minimization in temperature was performed on the expression Sj| log (pi — m log T — const |, 
where (pi is the set of 1000 solar model fluxes. Absolute values were considered to reduce slightly 
the weight of outliers. The inferred temperature exponents are not sensitive to the precise way 
the minimization is achieved. The best-estimate exponents given in Table are the average of the 
exponents computed by minimizing the residuals in either flux or temperature. The uncertainties 
presented span the range of the two solutions. 
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The last column of Table || gives the temperature exponents obtained in this section for 
a characteristic one-zone central temperature of Tc = 14 x lO^K. The derived exponents are 
not particularly sensitive to the assumed characteristic central temperature, Tc. Figure |^ 
shows the dependence of the exponents obtained with the one zone model for Tc in the range 
12 X lO^K to 16 X lO^K. 



A. pp and pep Temperature Exponents 

In the one-zone approximation, the measured solar luminosity can be written 

where V is volume, R is the rate of a nuclear fusion reaction per unit volume, R33 corresponds 
to the nuclear reaction ^He + '^He — > "^He + 2p (reaction 4 of Table D, -R34 corresponds to 
the reaction ^He + ^He — > ''Be + 7 (reaction 5 of Table |), and e is the amount of energy 
released by the fusion cycles corrected for neutrino losses. The rates of both R33 and R^^ 



increase rapidly with temperature; -R34 increases somewhat faster ||T8[. For the illustrative 
purposes of the one zone model, the three percent difference between £33 and €34 has been 
neglected in writing Eq. (H). We have also neglected the small contributions (less than a few 
percent in standard solar models) to the total luminosity of the fusion reactions associated 
with the CNO and ®B neutrinos. 

The flux of pp neutrinos at earth is 

'Pipp) = -^i2R33 + R34). (3) 

The factor of two appears in Eq. (^ because two pp neutrinos are produced in the first 
branch of the chain; a pp reaction is required to make each of the two ^He nuclei. In the 
second branch of the pp chain, the alpha particle acts as a catalyst and therefore only one 
pp reaction is needed. 

Substituting Eq. (□) into Eq. (H), one has 
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1 r?/" 1 

Since i?34 increases with temperature, the temperature dependence of the pp neutrino flux 



is negative, as was found first in the detailed evolutionary solar model calculations |T7 
The appropriate functional relation is therefore 

cj^ipp) = h-aTf (5) 

Here h = 6.5 x 10^''cm~^s~^ and 

, ciln[n(3He)n(4He)(3,4)] 

" 'ff' = din in ■ 

We make use of the convenient notation in which the reaction cross section times velocity 
averaged over a Maxwell-Boltzman distribution is represented, for the interacting nuclei i 
and j, by pointed brackets, (see Eq. 3.9 of reference [0). Thus the rate of the 

^He + "^He reaction is represented by i?34 = ?T,(^He) ?T,(^He)(3, 4). 

Both, the equilibrium number density of n{^B.e) and (3,4) are strong functions of tem- 
perature. The density can be found ||18|,^ by solving the equilibrium rate equations. At 
temperatures representative of the interior of the Sun, a good approximation for the number 
density of '^He is 



nCRe) ^ n(H) 



(U> p) 



The rate of the pp reaction is written Ri i = ra(-'^H)^(l, l)/2, where oc T~^/^ exp(— Tj ,,), 
Tij = 3Eo^ij/kT, and Eq is the most probable energy of interaction [l^. We note that r is 
proportional to T^^^'^. 

The value of m' can be calculated from Eq. (||) . 

Using this one-zone model, we can estimate the value of m' as well as the value of m, 
which is traditionally used to describe the temperature dependence [see Eqs. (P and (|])]: 

m'{pp) ^ 11; m{pp) ^ - m'{pp){aT"'' /b) ^ -OMm'ipp) ^ -0.9, (8) 
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where we have evaluated the exponents at the one-zone characteristic central temperature, 
= 14 X lO'^K. We have also taken aT"^' /h ^ to be equal to 0.08, as predicted 

by detailed numerical solar models or, less precisely, by the one-zone model. The preferred 
form for the temperature dependence of the pp neutrino flux is therefore 

(j){pp) oc [1-0.08(T,/T,,ssm)'"'], (9) 

where Tc,ssm = 15.64 x lO^K. 

The scaling exponent for the pp neutrino flux that is derived from the one-zone model 
agrees reasonably well with the exponent obtained from precise solar models (see Table ||) . 
The Monte Carlo study of 1000 detailed solar models yields m' = 13.0 (m = —1.1). The one- 
zone model yields m' = ll(m = —0.9). The one-zone model underestimates the exponent 
of the temperature dependence of the pp neutrino flux by about 20%. 

The exponent for the pep neutrino flux can be obtained with the aid of the analysis of 
the pp neutrino flux given above, since the rate of the pep reaction (reaction 2 of Table |I|) 
depends upon the rate of the pp reaction as R{pep) oc T~^^'^R{pp). Therefore, we can 
write 

<i>{pep) oc T-V2(T-) = T-^-\ (10) 

In the one-zone model, the pep neutrino flux, like the pp neutrino flux, scales like a negative 
power of the central solar temperature, as is found in the detailed solar model results. The 
numerical scaling derived from the detailed solar models has a rather large uncertainty, 
m{pep) = —2.4 ± 0.9, cf. Table Moreover, the power-law exponent is larger for the pep 
neutrino flux than it is for the pp neutrino flux, which follows from the fact that the pep rate 
divided by the pp rate is proportional to the modulus squared of the electron wave function 
near the two protons. The probability density of the electron is inversely proportional to 
the electron velocity I^Sf, which is itself approximately proportional to T^/^. 



Although the pep flux could be written in a form similar to Eq. (|^) for the pp neutrinos, 
we have chosen for simplicity to represent this minor component of the solar neutrino flux 
as a single power of Tc. 



Figure |^ shows that the derived scahng exponents for the pp and pep neutrino fluxes do 
not depend strongly on the value of the temperature, T^, which is assumed to characterize 
the solar interior in the one-zone model. 

We conclude this subsection by summarizing the main physical insight. The often-quoted 
dependence of the pp neutrino flux on a negative power of the temperature {(f){pp) oc T~^) 
results from the fact that as the central temperature gets larger, an increasing number 
of the completions of the pp chain proceed through the ^He + ^He reaction (reaction 5 of 
Table |). A complete fusion reaction of four protons being converted to an alpha particle 
via the '^He + ^He reaction involves the production of only one pp neutrino (cf. Table 1). On 
the other hand, a fusion of four protons via the ^He -|- ^He reaction (reaction 4 of Table |ID 
produces two pp neutrinos. The ^He -|- ^He reaction predominates at lower temperatures. 
Thus the pp flux is larger at lower central temperatures. 



B. ^Be and ^B Temperature Exponents 

The dependences of the other neutrino fluxes upon the central solar temperature can all 
be derived as simple power laws, oc T™. The power law exponent, m, can be obtained 
from the one-zone model. 

The flux of ^Be neutrinos can be calculated from the rate equation 

RCBe + e-) oc n(e)n(^Be)(e-, ^Be), (11) 

where R(^Be + e~) is the rate at which ''Be captures electrons in the solar interior (i.e., 
reaction 6 of Table |I|). Here (e~,'^Be) is, as before, the reaction cross section times velocity 
averaged over a Maxwell-Boltzman temperature distribution. The ^Be electron-capture 



reaction is much faster |T8[ under solar interior conditions than the competing proton-capture 
reaction (reaction 8 of Table |). The electron-capture rate is, in equilibrium, essentially 
equal to the rate of production of ^Be. Moreover, the production rate of ^Be is the same as 
the reaction rate i?34 that determines the temperature exponent m' of the pp reaction [see 
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Eq. (^) and the following discussion]. Therefore, the ^Be electron-capture reaction has the 
same scaling index, m(^Be) = 11, that was derived in Eq. for the pp reaction. Thus 

cj>CBe) oc Tl\ (12) 

which is approximately 10% larger than the value of m(^Be) = 10 that is obtained from the 
detailed numerical models (see Table 

The temperature dependence of the ^B neutrino flux can be calculated in a similar 
manner. The ^B neutrino flux results from a rare branch of the pp chain (reaction 8 of 
Table |) in which ^Be captures a proton rather than an electron. Therefore the '^B neutrino 
flux can be written as 

0(«B) cx ^rBe + e-)^|^. (13) 

Substituting = 14 x lO'^K in Eq. (|TD, we flnd: 

0(«B) oc Tf . (14) 

The temperature dependence is much stronger for ®B neutrinos than for ''Be neutrinos 
because the electron capture rate depends only weakly on temperature (essentially like 
T^^/^ @)' proton capture rate increases rapidly with temperature (like all strong 

interaction fusion rates). 

The derived scaling exponent for ^B, m(^'B) = 25, is fortuitously close to the value of 
m = 24 that is obtained by fltting to the detailed numerical models (cf. Table |T|). 

Figure |^ shows that the derived scaling exponents for the ^Be and *B neutrino fluxes 
vary by less than ±10% as the assumed characteristic central temperature Tc varies between 
12 X lO^K to 16 X lO^K. 



C. CNO Temperature Exponents 

The temperature dependence of the CNO neutrino fluxes can also be estimated simply. 
For the major part of the CNO cycle which leads from ^^C to ^^N, the slowest reaction is 
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"^^N(p, 7)^^0 [p^p4| . Therefore, the neutrino fluxes that are produced in this part of the 
cycle, from ^^N and ^^O, will both have approximately the same temperature dependence. 
(Slight differences will occur due, e.g., to non-equilibrium effects not accounted for in our 
static one- zone model.) The total number of CNO atoms is approximately constant and 
mostly in the form of ^^N. 

The temperature exponents for the ^^N and ^^O neutrino fluxes can be derived from the 
rate of the ^^N(p,7)^^0 reaction by calculating the logarithmic derivative of the reaction 
rate with respect to temperature. Thus m= [r("'^^N(p, 7)^^0) — 2]/3. The temperature 
dependence of the and ^^O neutrino fluxes is therefore: 

<PC^N), <PC'>0) oc Tf . (15) 

The exponents derived above are in reasonable agreement with the scaling laws obtained 
from the (non-equilibrium) detailed solar models, which are m(^^N) = 24 and m(^^O) = 27 
(see Table 0). 

Finally, we calculate the temperature dependence of the rare ^^F neutrino flux. The 
slowest reaction involved in producing the ^^F neutrinos is ^^0(p,7)^^F . The temperature 
dependence of the ^^F neutrinos can be calculated by analogy with the calculation for the 
^^N and ^^O neutrinos. In the derivation for the ^"^F neutrinos, we consider the ^^0(p,7)^''F 
reaction instead of the ^^N(p,7)^^0 reaction. The scaling law derived in this way is 

<j^C'F) oc Tf . (16) 

This result is in satisfactory agreement with the exponent obtained from the 1000 solar 
models, which is m(^^F) = 28 (see Table HI). 

Figure ^ shows that the derived scaling exponents for the CNO neutrino fluxes vary 
by less than ±10% as the assumed characteristic central temperature varies between 
12 X lO^K to 16 X lO^K. 
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IV. SUM RULE FOR TEMPERATURE EXPONENTS 



The luminosity of the sun can be expressed in terms of the individual neutrino fluxes, 
0i, as follows [0,0,13,11: 

Lq oc T,iei4>i, (17) 

where for neutrinos from the pp chain 

ei = 13.366 MeV- (gi). (18) 

Here (g,) is the average energy loss to the star from neutrinos of type i; the values of can 
be obtained from Table 3.2 of Ref. |T^. Equations (|T^ and ( [IS| ) assume that all the nuclear 



fusion reactions are in equilibrium, which is a reasonably accurate approximation for all but 
the CNO neutrinos and, in the outer region of the solar core, the reactions that produce and 
destroy ^He. For the ^^N and ^^O neutrinos, Eq. ([TsD does not apply and the values must be 
calculated separately from Table 3.3 of Ref. |jT8|. For these neutrinos, e(^^N) = 3.457 MeV 
and e(^^O) = 21.572 MeV. (Equation (18) would also apply to CNO neutrinos if they were 
in complete equilbirium. ) 

Fogli and Lisi [|14| pointed out that Eq. (0), when combined with the fact that the 



present-day solar luminosity is a known constant, implies a sum rule on the temperature 
exponents. In our notation, the Fogli-Lisi sum rule is 

Ejejmi^j = 0. (19) 

The Monte Carlo exponents given in Table H satisfy the Fogli-Lisi sum rule to an ac- 
curacy of 5% or better (i.e., (Sjejmi0j)/(Ejei0j) is less than 5%) . For the one-zone model, 
the sum rule is only satisfied to an accuracy of ~ 20%. (If only pp and ''Be neutrinos are 
considered, the sum rule is satisfied by construction in the one- zone model to the accuracy 
of our numerical approximations.) Violations of the sum rule in the one-zone model are 
caused primarily by the fact that different neutrino fluxes are produced in different temper- 
ature regions of the sun. If high precision is required, the solar neutrino fluxes cannot be 
parameterized by a single temperature. 
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V. CORRELATIONS BETWEEN NEUTRINO FLUXES 



In the previous sections, we have concentrated (as have most other investigations of 
this subject) on the average dependence of individual neutrino fluxes on the central solar 
temperature. In this section, we focus on the correlations that occur between the deviations 
of different neutrino fluxes from their average values. Previously we asked: On average, 
how strongly does a particular neutrino flux depend upon temperature? In this section, we 
ask: If one neutrino flux is larger than its average value by a specified amount, is a second 
neutrino flux larger or smaller than its average and, if so, by how much? 

The one-zone model predicts the relative magnitude and the relative phase of the frac- 
tional changes, A0/0, of the "^Be and the pp neutrino fluxes. Here 

A(j)i = (pi- 4>i,ssM, (20) 

where 0j,ssM is the standard value of the ith neutrino flux computed for the best input 
parameters and input physics. Since the temperature dependence of both the ^Be and the 
pp neutrino fluxes are, in the approximation in which we are working, governed by the 
rate, of the ^He + "^He reaction, the fractional changes in the fluxes are expected to 
be proportional to each other. The proportionality constant can be derived by comparing 
Eq. (D, Eq. (|), and Eq. (|Tg). We find 



A0(^Be) 
^(7Be) 



SSM 



(j){PP) SSM ^ ^ 



^Hpp) ^21) 



(j){pp) 



SSM 



'('^Be)ssM 

Figure ^ shows, as the one-zone model predicts, that the slope of the A(f)(^Be) / (f)(^Be) 
versus A0(pp)/0(pp) relation is negative and the magnitude of the slope is ~ —10. A closer 
study of Figure § reveals that the slope, a, obtained with the 1000 numerical models used 
in the Monte Carlo study is aMonte Cario ~ —9, whereas the slope predicted by the one-zone 
solar model is aone-zonc ~ —13. The difference between the slope obtained with the Monte 
Carlo study and the slope found with the one-zone model reflects the same imprecision in 
the one-zone model that was found earlier in Section [111 A| and Section [111B|. The slope that 
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is relevant for Figure ^ is the ratio of the ^Be and pp temperature exponents for neutrino 
fluxes; these exponents are predicted by the one- zone model to be (see Table |ID, respectively, 
10% too large and 20% too small relative to the detailed models. 

The one-zone model also predicts the average linear relation between the fractional 
changes, A0/0, of the ^B, ''Be, and pp neutrino fluxes. Figure ^ shows, in the top panel, the 
fractional changes in flux for ^B neutrinos versus the fractional changes in flux for the ''Be 
neutrinos, for the 1000 detailed solar models. The bottom panel of Figure |] shows fractional 
changes in flux for '^B neutrinos versus ''Be neutrinos. 

The relevant slope in a plot of A0j/0j versus A(f)j/(j)j is just the ratio of the corresponding 
temperature exponents for 0j and (pj that are given in Table ||. Therefore, the predicted 
one-zone model slope, a, for *B versus ^Be is Oone-zone ~ —2. 3, which is very close to the 
value of ttMonte Carlo ~ —2.4 fouud in the Monte Carlo study. The one- zone model predicts a 
somewhat too steep dependence of fractional changes in *B neutrino fluxes versus fractional 
changes in ''Be neutrino fluxes, namely, aone-zone ~ 28, versus aMonte Cario ~ 22. 

Figure ^ shows a large scatter in the relation between fractional changes of the ^B neutrino 
flux and either the ''Be or the pp neutrino flux. This lack of tightness in the relations shown in 
Figure ^ results ultimately from the fact that, in all modern solar models, the ^B-producing 
reaction, reaction 8 Table |, is rare and does not influence signiflcantly the structure of 
the sun. In fact, the largest uncertainty in the model calculations of the ^B neutrino flux 
is caused by the uncertainty in the experimental value for the low-energy nuclear cross 
section of reaction 8; the value of this cross section has essentially no effect (< 0.1%) on the 
calculated rates of the other nuclear fusion reactions. 

VI. CONCLUSION AND DISCUSSION 

The most interesting result of our study is the understanding it provides of the negative 
temperature dependence of the pp neutrino flux. The empirical fact that the pp flux decreases 
with increasing central temperature, contrary to the trend found with all other solar neutrino 
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fluxes, has been known since 1988 |1T7|JT8[] , but has not previously been explained physically. 



At first glance, this negative temperature dependence is counter-intuitive. 

In Section [111 A|, we show that the negative temperature dependence is a simple con- 



sequence of the fact that at higher temperatures only one pp neutrino is produced per 
(approximately) 25 MeV communicated to the star (via fusion of four protons), whereas at 
lower temperatures two pp neutrinos are produced per 25 MeV. In other words, at lower 
temperatures the ^He-^He fusion termination reaction (which requires two pp reactions) pre- 
dominates whereas at higher temperatures the ^He-^He reaction is faster (and requires only 
one pp reaction). The total energy per unit time communicated to the star must equal 
the observed solar luminosity, independent of the assumed central temperature. Thus, as 
the temperature increases and more of the nuclear fusion is accomplished by the ^He-^He 
reaction, fewer pp neutrinos are produced (and more ^Be and ®B neutrinos are created). 

In order to obtain a simple physical understanding of the temperature scalings and the 
correlations between the different neutrino fluxes, we have adopted a one-zone model for 
the interior of the sun. This model is characterized by a flxed central temperature, T^, 
and a total luminosity that is equal to the observed solar luminosity. Given the emphasis 
in the current literature on calculating ever more precise solar models, with hundreds of 
different mass shells, it is gratifying and surprising that the one-zone model accounts semi- 
quant it at ively for some of the most often used results of the detailed model calculations. 
Moreover, the one-zone model predicts the average correlations found between the different 
neutrino fluxes, a bonus in insight that was not possible to anticipate without detailed study 
of the simple model. 

Figure |I] and Figure |^ display the dependence upon central temperature of the 1000 
detailed solar models used in the Bahcall-Ulrich Monte Carlo study ||T7| , p!8| of theoretical 



uncertainties in the predicted solar neutrino fluxes. We have used these data to determine 
average temperature exponents, m, for all of the solar neutrino fluxes, where by assumption 
oc T"^. The exponents determined here are obtained by a formal best-fltting technique 



and are to be preferred to the previously-estimated exponents |117| , ^8| inferred less formally 
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from these same data; the previously-estimated exponents have been widely used in the 
literature. We have also estimated, for the first time so far as we know, uncertainties in the 
inferred temperature exponents. 

Figure ^ and Figure ^ show the correlations, found in the Monte Carlo study, between 
the different individual neutrino fluxes. These correlations reflect the fact that when one 
neutrino flux is increased or decreased, there is likely to be a corresponding change in the 
values of the other neutrino fluxes. These correlations must be taken into account when 
comparing the results of theoretical solar model calculations, including their uncertainties, 
with solar neutrino experiments. The only precise way to include the correlations displayed 
in Figure ^ and Figure || is to use the complete set of calculated neutrino fluxes in the 
theoretical analysis(cf. p^). Various practical approximations to this rather cumbersome 
method have been discussed in the literature (see, for example, |]6|J^JT^). 

The temperature exponents calculated with the aid of the one-zone model agree with 
the exponents inferred from the Monte Carlo study of precise solar models to an accuracy 
of 20% or better for the three most important solar neutrino fluxes: pp, ^Be, and ^B. The 
results are shown in Table which compares the exponents calculated with the one-zone 
model with the results obtained from the detailed solar models. Figure ^ shows that the 
scaling exponents calculated in the one-zone solar model are not strongly dependent upon 
the assumed characteristic central temperature, (taken here to be = 14 x lO^K). 

The quantitative agreement between the results of the one-zone model and the detailed 
models is impressive given the fact that the temperature exponents vary from m ~ — 1 for 
pp neutrinos to m ~ +24 for *B neutrinos. 

The physical insight provided by the one-zone model suggests a new form for the tem- 
perature dependence of the pp neutrino flux, which is given in Eq. (^). In this form, the 
variation of the pp neutrino flux is, for all temperatures, consistent with the observed solar 
luminosity, since it was derived by considering the relation between the solar luminosity, 
Eq. (2), and the pp neutrino flux, Eq. (3). Moreover, the formula for the pp neutrino 
flux, Eq. (^, provided by the one-zone model makes explicit the close correlation between 
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the ^Be and pp neutrino fluxes that is manifest in Figure |[ The expression used in this 
paper for the pp neutrino flux, Eq. (|^), was derived by considering the relation between the 
solar luminosity, Eq. (^), and the pp rate, Eq. (^. Physically, the strong correlation exists 
because the ^Be neutrino flux is proportional to the rate of reaction 5 of Table | and the pp 
neutrino flux is proportional to a constant minus the rate of reaction 5 (if we neglect the 
small contribution from CNO neutrinos). 

The one-zone model also accounts quantitatively for the average correlation, shown in 
Figure |^, between the ^B and ''Be neutrino fluxes, and between the ^B and pp neutrino 
fluxes. 

No simple model can, however, account in detail for the scatter in the correlation plots 
shown in Figure ^ and Figure ^. The Monte Carlo experiments simulate uncertainties 
in many different parameters; the power-law flts in the flgures represent only the average 
response of the neutrino fluxes to the changes in all the individual parameters. For analyses 
requiring a precise assessment of the correlations between the different neutrino fluxes, a 
Monte Carlo study of detailed solar models is required. 

What have we learned from this study? Improved temperature exponents for the neutrino 
fluxes are now available, with estimates for the uncertainties in the exponents. A static one- 
zone model of the sun accounts for the essential features of the temperature scaling of the 
neutrino fluxes and even describes well the average correlations between the fluxes. The 
model does not provide a precise description of the temperature dependences nor of the 
correlations between the different fluxes. The exponents derived from the one-zone mode 
model do not satify precisely the sum rule derived from the measured solar luminosity. 

The fundamental reason that the one- zone model does not account accuractely for all of 
the known results is that in precise solar models each neutrino flux is produced in a different 
range of temperatures. One cannot represent the results of different temperature ranges by 
a single parameter, Tc. 

In the future, a new Monte Carlo study must be undertaken to determine the temperature 
scalings and the correlations between the neutrino fluxes when, as required by helioseismo- 
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logical measurements |3^, diffusion is taken into account in the solar model calculations. 
The analysis of Bludman suggests that the effects of diffusion may alter the inferred 
temperature exponents by a non- negligible amount when compared to the values given in 
this paper, which are obtained from detailed solar models that do not include diffusion 
r?| . A Monte Carlo study is now underway that will create 1000 solar models that include 



diffusion and other recent refinements of the stellar model 31 
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FIGURES 

FIG. 1. The pp and pep Neutrino Fluxes as a Function of Central Solar Temperature. The top 
panel shows the pp neutrino flux versus the central solar temperature (expressed in units of lO^K). 
The lower panel shows the pep neutrino fluxes versus central temperature. The circles correspond 
to 1000 representations of the standard solar model calculated in a precise Monte Carlo study of 
the uncertainties in the standard model solar neutrino fluxes [|l^ . The two plotted lines represent 
a range of acceptable fits to the numerical data, which correspond to the indicated power-law 
dependences upon central temperature. The functional form for the pp reaction, 1 — 0.08{T/T^g^) 
is discussed in the text. 

FIG. 2. The "^Be, ^B, ^^N, and ^^O Neutrino Fluxes as a Function of Central Solar Temperature. 
The four panels show how the four different neutrino fluxes depend upon central temperature. The 
circles correspond to 1000 numerical solar models that were computed with different input data 



|17|. The plotted lines represent a range of acceptable fits to the numerical data, which correspond 



to the indicated power-law dependences upon temperature. 

FIG. 3. The Power-Law Exponents. The figure shows how the calculated power-law 
exponents(m, where neutrino flux oc T^) depend upon the characteristic central temperature of 
the one-zone model. For example, m{pp) varies between —0.9 to —0.8 as the central temperature 
is varied between 12 x lO^K to 16 x lO^K. 

FIG. 4. Correlation of the "^Be and the pp neutrino fluxes. The figure shows the strong corre- 
lation, predicted by the one-zone model, between the fractional changes in the '^Be neutrino flux 
and the fractional changes in the pp neutrino flux as calculated in the 1000 numerical solar models 



in the Monte Carlo experiment of reference |17|. Here A(/>j = (pi — (pi^sSM 



FIG. 5. Correlations between the ^Be, ^B, and pp neutrino fluxes. The figure shows the moder- 
ate correlation that exists between the fractional changes in the ^B neutrino flux and the fractional 
changes of either the ''Be or the pp neutrino flux. Here Acpi = — 4>i ssm • 
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TABLES 



TABLE I. The Principal Reactions of the pp Chain 



Reaction 


Reaction 


Neutrino Energy 




Number 


(MeV) 


rt 4- rt > ^TT -1- -1- ij 
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J. 


f) D tn D 4 
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1.4 


2tT _(_ ^ 3-tTp _|_ ^ 
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^He + ^He "^He + 2p 


4 
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^He + ^He ^ ^Be + 7 


5 




then 






e" + '''Be '''Li + 


6 


0.86, 0.38 


^Li + p ^ ^He + -^He 


7 




or 






p + ^Be ^ ^B + 7 


8 




^B ^ ^Be + e+ + i/e 


9 


to 14 



TABLE II. Temperature exponents for solar neutrino fluxes. We recom- 
mend a new functional dependence for <i>{pp) oc 1 — 0.08(T/Tssm)"*' as 
discussed in the text. All other exponents are given for the functional form, 
oc T^. For the one-zone model, we assumed a characteristic temperature 
Tc = 14 X lO^K. 
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